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Abstract. We present a rigorous determination of the critical value of the 
ground-state quantum Ising model in a transverse field, on a class of planar 
graphs which we call star-like. These include the star graph, which is a junction 
of several copies of Z at a single point. Our approach is to use the graphical, 
or FK-, representation of the model, and the probabilistic and geometric tools 
associated with it. 



1. Introduction 

The Hamiltonian of the quantum Ising model with transverse field on a finite 
graph G — {V, E) is the operator 

(1) H=-\x y: -iM''-^T.-i'' 

e—xy^E 

on the Hilbcrt space H = ^^.^y C^. Here the Pauh spin-1/2 matrices 

and we use as basis for each copy of in H the vectors l-l-)^; — (a) and \—)x = ( i ) ; 
also, X,6>0 are the spin-coupling and external-field intensities, respectively. Let 
/3 > denote the inverse temperature, and define the positive temperature states 

(3) PgAQ) = y^trie-P"Q), 

where Zc{P) = tr{e^^^) and Q G C^^^. Also define the ground state to be the limit 
Pg of pg,/3 as /? — + oo. If Gn is an increasing sequence of graphs tending to an infinite 
graph S, then we may also speak of infinite- volume limits ps,f3 = lim„_>oo PG„.f3 and 
PS = lim„_>oo Pg„ • The existence of these limits is discussed in [2] . 

In this article we will use the fk- or random-cluster representation of the ground 
state, see for example [12] and references therein. Details will be provided in the 
next section, but roughly speaking the FK-representation may be considered as a 
limit of "discrete time" random-cluster models on 5 x (eZ) as e J, 0. This is related 
to the well-known mapping of the quantum Ising model onto the classical Ising 
model in one dimension higher [17], and the FK-representation of that model [8]. 
The relevance of this representation is that it relates the occurrence of long range 
order in the ground state to the existence of infinite percolation paths in 5 x R; 
here we say that the model exhibits long range order if for all x, the correlation 
function 

(4) G(x,y) = p5((Ti3Vf ) 
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is bounded below by a positive function of x. There is a critical value of the ratio 
\/5 above which the model exhibits long range order, and below which it does not. 

The main result of this article is a rigorous determination of the critical ratio for 
a certain class of planar graphs S (see Definition 2). This extends the calculation 
for the graph S' = Z, to, amongst other graphs, the star graphs which is the junction 
of several copies of Z at a single point. See Figure 1. A special case of our main 
result (Theorem 8) is therefore the following. 

Theorem 1. The critical ratio for the ground state quantum Ising model on the 
star graph is A/5 = 2. 



Figure 1. The star graph has a central vertex of degree /c > 3 
and k infinite arms, on which each vertex has degree 2. In this 
illustration, k = A. 

In other words the critical ratio is the same for the star as for Z; this is to be 
expected since the star is only locally diflferent from Z. We emphasise, however, 
that the class of graphs for which we prove this result contains many more graphs 
than just the star. 

The quantum Ising model on Z has been thoroughly studied, and the critical ratio 
\/5 = 2 has been computed for this model in for example [16]. See also [17] and 
references therein. These calculations have relied on matrix methods and techniques 
such as Jordan-Wigner transformation. Recently, in [5], sharpness of the phase 
transition, and hence exponential decay of correlations below the critical point, 
was established rigorously for G = with any d > 1, using graphical methods 
similar to the corresponding proof [1] for the classical Ising model. Combining this 
result with duality arguments analogous to the classical two-dimensional random- 
cluster model [8], this gives another proof that the critical ratio is A/J = 2 for this 
model, using only tools from stachastic geometry (see [5] for details). One aim of 
this paper is to extend and illustrate the graphical methods, and show how they 
can be applied to a wider range of structures than just Z. The Ising model on 
the star-graph has also recently arisen in the study of boundary eff'ects in the two- 
dimensional classical Ising model, see for example [14, 15]. Similar geometries have 
also arisen in different problems in quantum theory, such as transport properties of 
quantum wire systems, see [6, 11, 13]. 



In this article we will let G = (V, E) be a star-like graph: 

Definition 2. A star-like graph is a countably infinite connected planar graph, 
in which all vertices have finite degree and only finitely many vertices have degree 
larger than two. 




2. Background and notation 
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Figure 2. A star-like graph G (left) and its line-hypergraph H 
(right). Any vertex of degree > 3 in G is associated with a "polyg- 
onal" hyperedge in H. 



Such a graph is illustrated in Figure 2; note that the graph of Theorem 1 is 
an example in which exactly one vertex has degree at least three. Fix a planar 
embedding G of G, and denote X = G x R; also let X = G x R := (V x R, £; x M). 
We will sometimes use X and X interchangably. Let O be a fixed but arbitrary 
vertex of G of degree two or more, which we think of as the origin. 

Recall that a hypergraph is a set W together with a collection F of subsets of W, 
called edges; a graph is a hypergraph in which all edges contain two elements. In 
our analysis we will use a suitably defined hypergraph "dual" of X: let H = {W, F) 
be the "line-hypergraph" of G, where W = E and the set {ei, . . . , e„} C E = W is 
in F if and only if 6i , . . . , Cn are all the edges adjacent to some particular vertex of 
G. Note that only finitely many edges of H have size larger than two. There is a 
natural planar embedding of H defined via the embedding G, in which an edge of 
size more than two is represented as a polygon. See Figure 2. Let Y = H x R and 
Y = H X M. 

Our configuration space fl will be the set of pairs cu = {B, D) where B C E x 'R 
and D <Z V X R are locally finite, which is to say that B fl ({e} x [—n,n]) and 
Dn{{v} X [-n, n]) are finite sets for slle e E,v eV and n e N. We think of B as 
a set of bridges and D as a set of deaths or cuts. There is a natural embedding of 
any oj € fl into X, where deaths are represented as missing points and bridges as 
"horizontal" lines connecting two "vertical" lines. See Figure 3 for an illustration 
of this when G = Z. Often we will identify w € with its embedding. Denote by 
d{-, •) the graph distance in G, and let A„ denote the set of points {v,t) and {e,t) 
where v gV has d{v, O) < n, e G E has at least one endpoint at distance at most 
n from O, and |t| < n. For each u € fl, we will employ two restricted embeddings 
Lo\ and w^, one "wired" and one "free". The free embedding a;° is simply the 
intersection of (the natural embedding of) w with A„. The wired embedding is 
defined by 

(5) = w° U {{v, t) : d{v, 0)=n + l, \t\ < n} U {(e, ±n) : e e G„}, 
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Figure 3. A configuration w on Z x M. Bridges are represented 
as horizontal line segments, and deaths as crosses. 



where we have taken the liberty to identify v gV and e £ E with their embeddings 
in G. In words, u)^ is obtained by tying together the top and bottom of a;°, as 
well as all bridges portruding from its "sides". We let the functions k^,kli_ : O — > N 
count the number of connected components of and w^, respectively. 

Equip Q, with the Skorokhod topology and the associated tr-algebra; the details 
of their definitions arc not immediately important, but may be found in [3] or [4]. 
Fix A, (5 > and let = fix^s be the probability measure on ft governed by a 
collection of independent Poisson processes Be on {e} x M, for e & E, and Dy on 
{v} X R, for V G V. Here each Be has intensity A, each D„ has intensity 6, and 
B = UeeB-Be, D = Uy^yDy. This fi is the space-time (or "continuum") percolation 
measure of [9]. We may now define the random-cluster probability measures. 

Definition 3. The random-cluster measure on A„ with parameters X,S,q > 
and boundary condition b G {0, 1} is the probability measure on O given by 



(6) 



Let 



(7) 9^ = $^((0, 0) lies in an unbounded component). 

The following basic facts may be proved in a conventional manner, as in [8, Theo- 
rem 5.5]; details for this particular model may be found in [4]. 

Proposition 4. Let ^ > 1. The weak limits := lim„_»oo exist, and enjoy a 

phase transition in the sense that there is pc = Pc{q) G (0, oo), depending only on 
q (and G), such that 0'' = if X/S < pc and 6*^ > if A/(5 > pc- We call pc the 
critical value of the random-cluster model on G x ffi. 

The relevance of the space-time random-cluster measures to the quantum Ising 

(or more generally quantum Potts) model is explained in [2] ; in particular the ground 
state quantum Potts model on G exhibits long-range-order iff the corresponding 
random- cluster model has 6^ > 0. Hence, to investigate the phase-diagram of the 
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quantum Ising model we will set q = 2 and focus on finding the critical value pc 
above which percolation occurs. 

Let us say a few more words about the "dual" Y of X. Given any configuration 
w G ri, one may associate with it a dual configuration on Y by placing a death 
wherever u has a bridge, and a (hyper)bridge wherever lo has a death. This is 
illustrated in Figure 4. More precisely, we let Oa be the set of pairs of locally finite 



Figure 4. Part of a configuration ui (solid) and its dual Wd (dashed 
with grey crosses) in the special case when G = Z. 

subsets of F X M and W xR, and for each w = {B, D) G Clwe define its dual to be 

LOd := (-D, B). As before, we may identify with its embedding in Y, noting that 
some bridges may be embedded as polygons. We let and \E''' denote the laws of 
Wd under and respectively. 

The case when G = Z is particularly important, and for this case we use the 
lower case symbols (p and ip in place of <1? and \E', respectively. When G = Z, the 
dual space Y is isomorphic to X, and we have the following result. Again the proof 
is similar to that for the discrete random-cluster model on Z^, but details for our 
model may be found in [4]. 

Lemma 5. If (j)'' have parameters q, A and 6, then the dual measures 

are random cluster measures with parameters q' = q, X' = qd and S' = X/q, and 

boundary condition 1 — 5. 

Recall that there is a partial order on fl given by {B' , D') — oj' > oj = {B, D) if 
B' D B and D' <Z D, and that an event A is called increasing if whenever G A 
and w' > oj then also uj' G A. Also recall that A is called a cylinder event if it 
only depends on a bounded region of X, which is to say that there is a bounded set 
A C X such that if a; = cj' on A then uj E A ii and only if co' G A. 

Definition 6. Let k be a probability measure on ^l. 

• We say that k is positively associated if for A, B any increasing cylinder 
events, k{A n B) > k{A)k{B). 

• Another probability measure ki on fl stochastically dominates k if for all 
increasing cylinder events A, we have ki{A) > k{A). We write ki > k. 
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• We say that k has the positivity property if for all £ > there exists a 
constant < c = c{e) < 1 such that for all e G E,v GV,t gB., 

(8) c < «;(no bridges in {e} x [t, f + e]) < 1 — c 
and 

(9) c < K(no deaths in {v} x [t,t + e]) < 1 — c. 

Proposition 7. Let q > 1. The measures i>J^, ^f^, '1/'' (6 = 0, 1) are positively 
associated and have the positivity property. Moreover, $^ > and > 

The proof of this is similar to the discrete random-cluster model and is omitted; 
full details may be found in [4]. 

3. The critical value 

We assume henceforth that q = 2. It is known that, if G = Z, the critical value 
Pc(2) = 2. The following is the main result of this paper. 

Theorem 8. Let G be any star-like graph. Then the critical value on G x M is 
Pc(2) = 2. 

In other words, the critical value for any star-like graph is the same as for Z. 

Simpler arguments than those presented here can be used to establish the analogous 
result when g = 1, namely that /9c(l) = 1- Also, the same arguments can be used 
to calculate the critical value of the discrete graphs G x Z when q = 1,2. 

Here is a brief outline of the proof of Theorem 8. First we make the straight- 
forward observation that /9c(2) < 2. Second, we use exponential decay and the 
GHS inequality to establish the existence of certain infinite paths in the dual model 
when \/S < 2. Finally, we show how to put these paths together to form "blocking 
circuits" in Y, which prevent the existence of infinite paths in X when X/6 < 2. 
Parts of the argument are inspired by [7] . 

Lemma 9. For G any star-like graph, Pc(2) < 2. 

Proof. Any star-like graph G contains an isomorphic copy of Z as a subgraph. Let 
Z be such a subgraph; we may assume that O G Z. Also we let (j)^, 4>^ denote the 
random-cluster measures on Z x'K. For each n > 1, let G„ be the event that in 
A„ there are no bridges between Z x R and its complement. Clearly each G„ is a 
decreasing event. It follows from a standard property of random-cluster measures, 
sometimes called the DLR-property, that $^(- | G„) = (j^ni-)- The proof of this 
uses standard techniques [8]; details for this model may be found in [4]. If A is an 
increasing cylinder event, this means that 

(10) cli{A) = ^l{A\Cn)<^i{A), 

i.e. < for all n. Letting n ^ oo it follows that cj)^ < If \/5 > 2 then 
(t>''{{0,0) oo) > so then also 

(11) $*((O,0) ^ oo) > 0, 



which is to say that Pc{2) < 2. 



□ 
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3.1. Infinite paths in the half-plane. Let us now establish some facts about 
the random-cluster model on Z+ x R which will be useful later. Our notation is as 
follows: for n > 1, 

S„ = {(a, t) e Z X R : -n < a < n,\t\ < n} 

^"^^^ Sn{m, s) = Sn + (m, s) = {{a + m,t + s) e Z+ X R : (a, t) e 

For brevity write T„ = S'„(n,0); also let d denote the boundary, 

(13) dSn = {ia,t) e Z X R : a = ±n or t = ±n} 

and dSn{m, ,s) = dSn + (m, s). For 6 = 0, 1 and A one of 5„, T„, we let denote 
the q = 2 random-cluster measure on A with boundary condition b and parameters 
X,6. Note that 

(14) (/)'' = lim (j)^ , i/' = lim ^'s • 
We will also be using the limits 

(15) (j/" = lim 0^ , V-^ = lim ipr ■ 

n — yoo " n — ^oo " 

These are measures on configurations w on Z+ x R; but according to our definition 
they cannot be random-cluster measures since the regions r„ do not tend to the 
whole of Z X R. However, standard arguments let us deduce all the properties of 
(/)"', tjj-f that we need. In particular ijj^ and t/)™ are mutually dual (with the obvious 
interpretation of duality) and they enjoy the positive association and positivity 
properties of Definition 6. 
Let W be the "wedge" 

(16) W = {{a,t) &Z+ xR:0<t<a/2 + l}, 
and write for the origin (0,0). 

Lemma 10. Let X/S < 2. Then 

(17) (0 oo in W) > 0. 

Here is some intuition behind the proof of Lemma 10. The claim is well-known 
with tp° in place of ijj^ , by standard arguments using duality and exponential decay. 
However, is stochastically smaller than ■0", so we cannot deduce the result 
immediately. Instead we pass to the dual f/)™ and establish directly a lack of blocking 
paths. The problem is the presence of the infinite "wired side"; we get the required 
fast decay of two-point functions by using the following result of [10], adapted to 
our model. 

Proposition 11. Let X/S < 2. There is a > such that for all n, 

(18) (f>s^{0^ dSn) <er"''. 

In words, the two-point function decays exponentially also in finite volume. 
Higuchi [10] proves a more general result for the discrete Ising model, but at- 
tributes to Aizenman the simpler result for that model. It was pointed out to us 
by Grimmett (personal communication) that the original proof may be shortened 
by using the Lieb inequality in place of the GHS inequality, and we present the full 
proof for our model using the Lieb inequality here. Apart from the Lieb inequality, 
the proof uses another fact known for the g = 2 Ising case but not for the general 
case g > 1, namely exponential decay in the infinite volume subcritical Gibbs state. 
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Proof. Let S'„ ^ Sn denote the "tall" box 

(19) 5„ = {{a, t) e Z X m. : < a < n,\t\ < n + 1}. 

We will use a variant of the random-cluster measure on Sn which has non-constant 
intensities for bridges and deaths, and also a process of ghost-bonds. To this end 
we create a new site g, which we think of as a "point at infinity", and let S(-),j(-) : 
Z X M ^ R and A(-) : (Z -|- 1/2) x M ^ M be bounded, nonnegative and measurable 
functions. Given independent Poisson processes of bridges and deaths of rates A(-) 
and 5(-), respectively, and of links to g of rate 7(-), we may define random-cluster 
measures as in Definition 3, where now any components connected to g are to be 
counted as the same. 

The particular intensities we use are these. Fix n, and fix m > 0, which we think 
of as large. Let A(-), 6{-) and 7m(-) be given by 



6{a,t) 



5, if (a, t) G Sn 
0, otherwise, 

A, if (a, t) e Sn and {a+l,t) e Sn 



(20) ^(« + V2,i) - I otherwise, 

{A, if exactly one of (a, t) and (a -t- 1, t) is in Sn 
m, a {a,t) G Sn\Sn 
0, otherwise. 

In words, the intensities arc as usual "inside" Sn and in particular there is no 
external field in the interior; on the left and right sides of Sn, the external field 
simulates the wired boundary condition; and on top and bottom, the external field 
simulates an approximate wired boundary (as m oo). We introduce another 
parameter r G [0, 1], and let cp''^ „ denote the random-cluster measure on S„ with 

intensities A(-), S{-), rjm{-). Note that t^'^^ „ and (f/g agree on events defined on Sn, 
for any m. 

Let X denote Sn \ Sn together with the k^ft and right sides of Sn ■ By the Licb 
inequality, proved for the space-time Ising formulation of the present model in [5] 
(see also [4]), we have that 
(21) 

^ln,n{0 ^9)< e'' / dx 0^,„(O ^ x)^Ux ^g)< e'' [ dx ^^,„(0 ^ x), 
Jx Jx 

since X separates from g. Therefore, by stochastic domination by the infinite- 
volume measure, 

(22) 4>l^^{0 ^g)< e'' [ dx </)°(0 ^ x). 

Jx 

All the points x G X are at distance at least n from the origin. By exponential 

decay in the infinite volume, as proved in [5] using similar methods to the discrete 
case [1], there is an absolute constant a > such that 

(23) 4>ln,ni^ 5) < e«^|X|e-"" = e^'iSn + 2)6""". 

Now let C be the event that all of Sn \ Sn belongs to the connected component of g, 
which is to say that all points on Sn \ Sn are linked to g. Then by the DLR-property 
of random-cluster measures the conditional measure (f>m,n{' I ^) agrees with (f>s^{-) 
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on events defined on Sn- Therefore 

4„(0 ^ dSr.) = ^l^^iO ^ dSn I C) = 0^,„(O ^g\C) 

(24) ^1 (0 ^ a) e^* 

^ ^ <TlIhI^ (8n + 2)e-"". 

Since 4>m ti(^) — > 1 as m — > oo we conclude that 

(25) 4„ (0 ^ dSn) < e^\8n + 2)6""". 

Since each (0 dSn) < 1 it is a simple matter to tidy this up to get the result 
claimed. □ 

of Lemma 10. Let T = {{a,a/2 + 1) : a € Z+} be the "top" of the wedge W. We 
claim that 

(26) ^0'^((n,O) ^ T in W^) < oo. 

n>l 

Once this is proved, it follows from the Borel-Cantelli lemma that with probability 
one under t/)™, at most finitely many of the points (n, 0) are connected to T inside 
W. Hence under the dual measure ip^ there is an infinite path inside W with 
probability one, and by the positivity- and positive association properties it follows 
that 

(27) V^(0 ^ oo in W) > 0, 
as required. 

To prove the claim we note that, if n is larger than some constant, then the 
event "(n,0) ^ T in M^" implies the event "(n,0) ^ 9S'„/3(n, 0)" . The latter 
event, being increasing, is more likely under the measure (pg than under 

But by Proposition 11, 

(28) 4„/3(n,o)((«>0) ^ 95„/3(n,0)) = <f>\^jO ^ dS^/,) < e-«"/^ 

which is clearly summable. □ 

3.2. Proof of the main result. We prove one more lemma about the half-plane 

before going on to the main result. 

Lemma 12. Let X/6 < 2. There exists e > such that for each n, 

(29) ((0, 2n + 1) ^ (0, -2n - 1) off r„) > e. 

Proof. Let L„ = {(a, n) : a > 0)} be the horizontal line at height n, and let £ > 
be such that ^'■'^(0 oo in W) > \fe. We claim that 

(30) V-^ ((0, -2n - 1) ^ L2„+i off T„) > ^fe. 

Clearly is invariant under reflection in the x-axis, and standard arguments [8, 
Theorem 4.19] imply that it is also invariant under vertical translation. Thus once 
the claim is proved we get that 

(31) ((0, 2n + 1) ^ (0, -2n - 1) off r„) > 

> V-^((0, -2n - 1) ^ L2n+i off r„ 

and (0, 2n + 1) ^ L-in-\ off r„) > {^ef , 

as required. See Figure 5. 
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L-2n-l 



Figure 5. Construction of a "half-circuit" in Z+ x M. With prob- 
abihty one, any infinite path in the lower wedge must reach the 
line L2n+i, and similarly for any infinite path in the upside-down 
wedge. Any pair of such paths starting on the horizontal axis must 
cross. 

The claim follows if we prove that 

(32) V^(0 ^ oo in R) = 0, 
where R is the strip 

(33) R={{a,t) ■.a>0,-2n-l<t<2n+l}. 

However, (32) follows from the positivity property of Definition 6 and the Borcl- 
Cantelli lemma, since the event "no bridges between {k} x [— 2n — 1, 2n + 1] and 
{fc + 1} X [— 2n — 1, 2n + 1]" must happen for infinitely many k with ^/'•''-probability 
one. To see this we can compare with an independent percolation measure, as 
in the proof of Proposition 11. We have that < ji, where has parameters A, 5; 
under fi the events above are independent, so 

(34) (0 <-> oo in i?) < /i(0 <-> oo in i?) = 0. 

□ 

of Theorem 8. We may assume that G ^ Z, since the case G = Z is known. Let 
\/5 < 2, and recall that G consists of finitely many infinite "arms", where each 
vertex has degree two, together with a "central" collection of other vertices. On 
each of the arms, let us fix one arbitrary vertex (of degree two) and call it an exit 
point. Let U denote the set of exit points of G. 

Given an exit point u € U, call its two neighbours v and w; we may assume 
that they are labelled so that only v can reach the origin O without passing u. 
If the edge uv were removed from G, the resulting graph would consist of two 
components, where we denote by J„ the component containing w. Let 
denote the marginals of ^'^,i>^ on X„ := J„ x M; similarly let denote the 

marginals of the dual measures. Of course is isomorphic to the half-plane graph 
considered in the previous subsection. By positive association and the DLR-property 
of random-cluster measures, $° < 0^ so letting n — > oo also < (f)^. Passing 

to the dual, it follows that > tp^ . The (primal) edge uv is a vertex in the 
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line-hypergraph; denoting it still by uv we therefore have by Lemma 12 that there 
is an £ > such that for all n, 



Here Tn{u) denotes the copy of the box Tn contained in Xu- Letting A denote 
the intersection of the events above over all exit points u, and letting Ai = Ai{n) 
be the dual event Ai — {u^ : co G A}, it follows from positive association that 
^^{Ai) > , where k — \U\ is the number of exit points. Note that Ai is a 
decreasing event in the primal model. The intuition is that on Ai, no point in 
Tn{u) can reach oo without passing the line {u] x [— 2n — 1, 2n + 1], since there is 
a dual blocking path in X„. 

Next let / denote the (finite) subgraph of G spanned by the complement of all 
the J„ for u G U, and let A2 = A2{n) denote the event that for all vertices v G I, 
the intervals {v} x [2n + 1, 2n + 2] and {v} x [—2n — 1, —2n — 2] all contain at least 
one death and the endpoints of no bridges (in the primal model). By the positivity 
property, there is rj > independent of n such that $°(j42) > rj. So by positive 
association $"(^1 0^12) > V^'^ > 0- On the event Ai nA2, no point inside the union 
of / X [—71, n] with (JueuTn{u) can lie on an infinite path. See Figure 6. Taking the 



Figure 6. The dashed lines indicate dual paths that block any 
primal connection from the interior to 00. Note that this figure 
illustrates only the simplest case when G is a junction of lines at 
a single point. 

intersection of the Ai{n) n A2{n) over all n, it follows that 

(36) $° (there is no unbounded connected component) > rje . 

The event that there is no unbounded connected component is a tail event. All 
infinite- volume random-cluster measures are tail-trivial (see [8, Theorem 4.19] or [4]), 
so it follows, whenever X/5 < 2, that 



(35) 



■^'^{{uv, ~2n - 1) ^ {uv, 2n + 1) off T„(m) in Xu) > e. 





(37) 



$0 



(0 <A 00) = 1. 



In other words, Pc{2) > 2. Combined with the opposite bound in Lemma 9, this 
gives the result. □ 
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